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ABSTRACT 

Strong evidence for the presence of a warped Keplerian accretion disc in NGC 4258 
(M 106) has been inferred from the kinematics of water masers detected at sub-parsec 
scales. Assuming a power-law accretion disc and using constraints on the disc pa- 
rameters derived from observational data, we have analyzed the relativistic Bardeen- 
Petterson effect driven by a Kerr black hole as the potential physical mechanism 
responsible for the disc warping. We found that the Bardeen-Petterson radius is com- 
parable to or smaller than the inner radius of the maser disc (independent of the 
allowed value for the black hole spin parameter). Numerical simulations for a wide 
range of physical conditions have shown that the evolution of a misaligned disc due 
to the Bardeen-Petterson torques usually produces an inner flat disc and a warped 
transition region with a smooth gradient in the tilt and twist angles. Since this struc- 
ture is similar to that seen in NGC 4258, we propose that the Bardeen-Petterson effect 
may be responsible for the disc warping in this galaxy. We estimated the time-scale 
necessary for the disc inside of the Bardeen-Petterson radius to align with the black 
hole's equator, as a function of the black hole spin. Our results show that the Bardeen- 
Petterson effect can align the disc within a few billion years in the case of NGC 4258. 
Finally, we show that if the observed curvature of the outer anomalous arms in the 
galactic disc of NGC 4258 is associated with the precession of its radio jet/counterjet, 
then the Bardeen-Petterson effect can provide the required precession period. 

Key words: galaxies: active - galaxies: individual: (NGC 4258) - accretion, accretion 
discs - masers - black hole physics - relativity 



1 INTRODUCTION 

The barred spiral galaxy NGC 4258, loc ated at a distance 
of 7.2 ±0.3 Mpc l|Herrnstein et alj|l99g| h presents a well- 
studied ac tive nucle a r regi on, classified as a Seyfert 1.9, ac- 
cording to lHo et al.l (|l997l h with a warped sub-parsec Kep- 
lerian water m aser accretion disc surrounding a supermas - 
sive black hole (jMivoshi et al.lll995l : iHerrnstein et al1 !2005) . 
iNeufeld fc Malonevl il99^ ) modelled the maser region as a 
warped viscous accretion disc that is illuminated obliquely 
by a central X-ray source, so that the disc heating is domi- 
nated by X-ray photons. They also suggested that the edges 
of the water maser region are determined by the thermody- 
namical conditions of the disc (and this has been corrobo- 
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rated recently bv lHerrnstein et al.| [2005l . Low-energy X-rays 
from the galactic centre are absorbed by a hydrogen column 
density of approximately 10 23 cm -2 , providing indirect evi- 
dence that th e warped accretion di sc is the X-ray absorber 
in NGC 4258 (|Fruscione et al.ll2005h . 

High-resolution interferometric image s at radio wave- 
lengt hs show the presence of a compact jet (jHerrnstein et all 
1 1997T ) . and of two radio spots located at 840 pc south and 1.7 
kpc north from the nucleus and identifie d as bow shock co un- 
terparts of the jet-ambient interaction |Cecil et al.ll200(ih . 

The maser disc, which extends roug hly between 0.14 
and 0.28 pc in relation to the nucleus (|Herrnstein et all 
2005), is count errot at ing with respect to the galactic 
disc of the host galaxy, in which kilo-parsec anoma- 
lous arms have been de tected at the Ha transition 
l|Courtes fc Cruvellierl Il96l[ ) and later at radio to X-ray 
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A detaile d study of the kinemati cs of the maser spots 
conducted by iHerrnstein et all (|2005h showed a deviation 
from Keplerian motion in their projected rotation curve of 
about 9 km s _1 . This deviation was successfully modeled as 
a Keplerian warped accretion disc, with a radial gradient in 
its inclination of approximately 0.034 mas -1 . 

Although the existence of a warped accretion disc in 
NGC4258 seems to be the best explanation for the observa- 
tional data, the nature of the physical mechanism responsi- 
ble for producing such a structure has not been unambigu- 
ously established. 

Tidal torques in a massive binary system were proposed 
bv IPapaloizou. Terauem fc Linl (| 199Sl } as being responsible 
for the disc warping (in their model, the companion object 
has a mass comparable to that of the disc). 

In addition, irradiation-driven torques prod uced by 
non-a xisymmetric forces due to radiation pressure l)Pringlel 
1996) were also in voked to explain the observ e d be - 
havi our. As noted by Malonev. Begelman fc Pringlel (|l996l ) 



and iMalonev. Begelman fc Nowakl 1 19981 ), if the accretion 
disc o f NGC 4258 is radiatively inefficient (jLasota et al.l 
1 19961 ; [G ammie. Naravan fc Blandfordl Il999l ). this mecha- 
nism cannot produce the inferred warp in NGC 4258, since 
the disc becomes radiatively unstable far from its maser re- 
gion. If the disc is radiatively efficient but the shear act- 
ing vertically does not have the same value as the az- 
imuthal shear, the critical radius beyond which the ra- 
diative instability acts will increase. For a viscosity pa- 
rameter a < 0.2 (|Shakura fc Sunvaevl [19731 ). the masing 
disc of NGC 4258 is stable against the radiation instability 
l|Gammie. Naravan fc Blandfordlll999l ). 

In this work, we analyze an alternative mechanism that 
could warp the accretion disc of NGC 4258 : the Bardeen- 
Petterson effect |Bardeen fc Petterson! Il975h . predicted in 
the framework of general relativity when the spin axis of 
a Kerr black hole is inclined in relation to the angular mo- 
mentum vector of the disc. A brief description of the effect is 
presented in Section 2. In Section 3, we present our accretion 
disc model, introducing its basic parameters, as well as their 
values constrained by the observational data of NGC 4258. 
The Bardeen-Petterson radius and the alignment time-scale 
between the accretion disc and the spinning black hole, cal- 
culated in terms of the black hole spin and for different disc 
parameters, are presented in Section 4. The possible connec- 
tion between the radio jet of NGC 4258 and the anomalous 
arms in the galactic disc is also discussed in Section 4 in 
terms of jet precession due to the Bardeen-Petterson effect. 
General conclusions are summarized in Section 5. 



mainly studied analytically and numerically in the con- 
text of quas i-periodic X-ray brightness oscillations in X- 
ray binaries dStella fc Vietr|ll99a iMarkovic fc Lamblll9"98l ; 
lArmitage fc Nataraianlll999h . The precession angular veloc- 
ity Olt due to the Lense-Thirring effect is given by (e.g., 
Wilkins 1972): 
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where G is the gravitational constant, c is the speed of light, 
r is the distance from the black hole and Jbh is the angular 
momentum of a Kerr black hole with mass Mbh, which is 
defined as: 
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GM£ 
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where a, (|a*| ^ 1) is a dimensionless parameter correspond- 
ing to the ratio between the actual angular momentum of 
the black hole and its maximum possible value. 

The alignment between the angular momenta of the 
Kerr black hole and the accretion disc is forced by a 
combination between the Lense-Thirring effect and the 
internal viscosity of the ac cretion disc. This is known a s the 
Bardeen-Petterson effect (jBardeen fc Petterson! Il975h and 
it tends to affect the innermost part of the disc due to the 
short range of the Lense-Thirring effect, while the disc's 
outer part tends to remain in its original configuration. The 
transition radius between these two regimes is known as 
the Bardeen-Petterson radius 7?bp and its exact location 
depends m ainly on the physical prope r ties of the accre- 
tion disc dBardeen fc Petterson) 197a: Kumar fc Pringlel 
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Basically, the Bardeen-Petterson radius is determined by 
comparing the time-scale related to the Lense-Thirring 
effect to that of warp transmission through the disc, giving: 
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2 THE BARDEEN-PETTERSON EFFECT 

2.1 The Bardeen-Petterson radius 

Frame dragging produced by a Kerr black hole causes 
precession of a particle if its orbital plane is inclined 
in relation to the equatorial plane of the black hole. 
This effect, known as Lense-Thirring precession, has been 



depending on whether the warp propagation occurs diffu- 
sively (-Rbp ) or through waves (R^ P ). In equations (3) and 
(4), V2 represents the viscosity along the normal to the ac- 
cret ion disc, and c s is t he sou nd speed in the disc. 

IPapaloizou fc Linl i|l995h showed that the transition 
from the diffusive to wave-like regime occurs at a radius 
Rt ~ Hd/a, where iJd is the scale-height of the disc and 
a is the dimensionless vi scosity parameter introduced by 
IShakura fc Sunvaevl (|l973l ). 
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2.2 Alignment time-scale 



3.1 Parametric model for the accretion disc 



The time-scale for alignment between the angular momenta 
of the bl ack h o le an d the accretion disc was first esti- 
mated by iReesI i| 19781 ), assuming that each mass element 
accreted by the black hole carries an orbital angular momen- 
tum corresponding to t hat at the Bardeen-Petterson radius. 
IScheuer fc Feilerl (Il996l ) obtained an analytic solution to the 
equations that control the warp evolution in the case of a 
disc with consta nt surface density and calculat ed the align- 
ment time-scale. iNataraian fc Armitagel (1 1999T) generalized 
the results found by Scheuer fc Feilerl l| 19961 ) to a power-law 
viscosity. These studies suggest that the alignment time- 
scale can be estimated by: 



align 
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( dJBH_\ " 
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where tp is the angle between the black hole spin axis and 
the direction perpendicular to the outer part of the warped 
disc. The time derivative of Jbh has the form: 



We will consider an accretion disc with total mass Md- In the 
case of an axisymmetric disc, its total mass can be calculated 
through: 



M d (r <: Rout) = 2tt 



T,(r)rdr. 



(9) 



Here 7? ou t is the outer radius of the accretion disc, 
while Rms = A ms (a„)i? g is the radius of the marginally 
stable orbit (assumed to be the inner radius of the disc), 
R g — GMbh/c 2 is the gravitational radius, and A ms (a*) = 
3 + A 2 T y/{3-A 1 ){3 + A 1 + 2A 2 ), with Ai = 1 + (1 - 



2N1/3 



[0+.. 
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+ (l-q«) 1/3 ] and A 2 = 3o| + A\ 



(e.g., iBardeen. Press fc~ Teukolsky|[l972l) . In this case, the 
minus and plus signs correspond to prograde and retrograde 
motion, respectively. 

The scale-he ight of the disc Hd can be esti- 



mated following Sakimoto fc Coronitil (1 1 98 lh (see also 
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where Ld is the angular momentum density of the accretion 
disc, given by: 



L d (r) = E(r)fi d (r)r 2 
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where Q d and E are respectively the angular velocity and the 
mass surface density of the accretion disc. The last quantity 
is obtained by integrating over the scale-height of the disc: 



r-Hd/2 

E(r) = / p(r,z)dz, 

-H a /2 



(8) 



where p is the mass density of the disc and z is the perpen- 
dicular distance from the midplane of the disc. 

Note that the alignment torque depends on an integral 
that should be performed from the Bardeen-Petterson radius 
to the outer radius of the accretion disc. Although the last 
quantity is not known accurately in many cases (even in 
the case of NGC 4258), its precise value is not important for 
our calculations since the alignment torque decreases rapidly 
with distance for the realistic disc models considered in this 
work. 



3 ACCRETION DISC OF NGC 4258 

The examination of the viability of the Bardeen-Petterson 
effect requires some knowledge about the physical charac- 
teristics of the accretion disc (such as its viscosity, surface 
density and scale-height), as well as the mass and spin of 
the black hole. In this section, we will present our paramet- 
ric accretion disc model for NGC 4258, using the available 
observational data to constrain the values of the parameters. 



H ss = Cs 2 /(vrGE d ) 



(12) 



This formulation takes into account possible effects re- 
lated to self-gravity. In the case of NGC 4258, maser ob- 
servations have shown that this effect is negligible (e.g., 
iHerrnstein et al.ll2005T) . so that we will assume hereafter 
#d = H nBg . 

The sound speed c a is given by (e.g.. lAbramowicz et al.l 

ll98Sfr : 



c s (r) 



,dlnfi d (r) vi (r)n d (r) 



din r 
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where P is the politropic index of the gas, which we have 
taken to be equal to 5/3, and v\ is the viscosity along the 
disc: 



v\(r) 



M 



27rE d (r) 



a! In fid (r) 
dlnr 



Qd(Rm: 

fid(r) 



where M is the accretion rate onto the black hole, assumed 
to be constant along the disc. 

The viscosity acting in the vertical direction, v 2 , is 
written in terms of v\ as v 2 ~ f(ct )v\, where f(a) = 
2(1 + 7a 2 )/[a 2 (4 + a 2 )] (|Ogilvie! mm ) . If the warp prop- 
agates diffusively along the d isc in the linear regime an d 
a < 1, we have f(a) w l/2a 2 (|Papaloizou fc Pringle|[l983l ). 
Since f(a) is substantially greater than unity for the usual 
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Figure 1. Cumulative mass of the accretion disc of NGC4258, 
normalized by the black hole mass, as a function of radius for 
the upper limit of Eq. The continuous line refers to a mass sur- 
face density of the disc decreasing with the square of the dis- 
tance (s = —2), while the dashed and dashed-dotted lines cor- 
respond respectively to s = — 1 and 0. The vertical dashed lines 
mark the inner and outer radii of the maser disc. The horizon- 
tal dashed line shows the observed upper limit for M^/A^bh 
l|Moran. Greenhill fc Herrnsteinlll999h . 



values of a, warping modes very effectively propagate in- 
ward the disc, where they are dissipated later. Note that 
in the strong non-linear regime (warps with very large am- 
plitudes), additional dissipation cause d by fluid instabilitie s 
might reduce /(a), so that v^jvx ~ 1 (iGammie et al.| [2000). 
Both possibilities were studied bv lLodato fc Pringld T2006h 
in their numerical simulations of the Bardeen-Petterson ef- 
fect. In this work, we are assuming that the linear regime is 
applicable for the accretion disc of NGC 4258. 



3.2 Constraining accretion disc model parameters 
for NGC 4258 from maser observations 

3.2.1 Mass of the black hole and accretion disc rotation 
law 

The black hole mass in NGC 4258 has be en determined 
from interferometric maser observations (Mivoshi et al 



19951: Im oran. Greenhill fc Herrnsteinll 999: Hcrr nstein et al 



Herrnstein et al.1 (|2005T ) 



The rotation curve of the mas er spots is compat- 
ible with Keplerian motions (e.g., iMivoshi et all 1 19951 : 



iHerrnstein et al]|2005h . Therefore, we will assume a (rela- 
tivistic) Keplerian angular velocity £Ik for the accretion disc 
of the form: 



n K (r) = 



GMbh 



3/2 



(15) 



Table 1. Accretion disc parameters for NGC 4258. 
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Upper limit. 

Highest value among all black hole spins (eM ^ 1.77 X 10~ 6 Mq yr _1 ). 
Lowest value among all black hole spins (a $J 0.2). 

Highest value among all black hole spin values (ibol/^Edd 1? 2.10 X 10 — 5 ). 



3.2.2 Surface density of the accretion disc 

We have adopted a power-law rad ial distribution for the 
accretion disc surface density (e.g., Papaloizou fc Terqueml 
1 19951 : lLarwoodlH997l : iNelson fc Papaloizou! l200Ch . 



(16) 



where Eo and s are constants to be determined from obser- 
vations. 

From the relation between mass and surface densities of 
the accretion disc (equation 8) , Eo can be expressed through: 



(17) 



where p = pman, p is the mean molecular weight, and ran 
is the mass of the hydrogen atom. 

The maser amplification is favoured w hen the part i- 
cle density n ranges from 10 8 to 10 10 cm" 3 (|Elitzurlll993 l. 
Thus, for NGC 4258, this condition must be fulfilled. A sta- 
tistical analysis of the vertical structure in the systemic 
masers of NGC 4258 h as suggested that H d /r < 0.002 
l|Herrnstein et al1l2005h . Taking into account this observa- 
tional constraint, equation (17) provides an upper limit on 
E , such that E;T x (s) < 1.85 X 10~ 4 (r = 0.14 pc/flg) 1 " 8 g 
cm -2 , calculated at the inner radius of the maser disc for a 
particle density of 10 10 cm -3 . 

On the other hand, using equation (11) for the scale- 
height of the disc, we obtain: 



20051): we will adopt M B h = 3.78 x 10 7 M Q , as inferred by s o 



2TvR e 



-0+1/2) 



(18) 



Using again the observational constraint Hd/r < 0.002 
at r = 0.14 pc in equation (18), we obtain a lower limit for 
E , denoted by Eg 1 ' 11 . 

Note that this equation was obtained under the assump- 
tion that viscosity is responsible for disc heating. However, 
the warped maser region is also heated by the X-ray con - 
tinuum from the central source l|Neufeld fc Malonevlfl995l ). 
which increases the magnitude of the local sound speed rel- 
ative to that calculated from equation (13). Consequently, 
it could change the value of E™ 111 calculated via equation 
(18). In this paper we will neglect the X-ray heating since 
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Figure 2. Accretion disc parameters of NGC 4258. Upper panels: The surface density at r = R g for three different power-law dependencies 
of the disc surface density (s = —2, —1 and 0) and the accretion rate multiplied by the accretion efficiency as a function of the black hole 
spin. The gray areas show the range of parameters that obey the observational constraints. Lower panels: The alpha parameter and the 
bolometric luminosity in units of the Eddington limit as a function of the black hole spin (same nomenclature as in the previous panels 
is used). 



its inclusion does not change substantially either our es- 
timates of the Bardeen-Petterson radius or the associated 
ali gnment timescale. The us e of an isothermal model as 
in iNeufeld fc Malonevl fl995) would change the Bardeen- 
Petterson radius by a factor in a range of about 0.2 - 2, 
while for the alignment timescale the correction would not 
be larger than a factor of two. 

From the upper limit of 10 6 Mq for the 
disc mass inferred from t he m as er observations 
llMoran. Greenhill fc Herrnsteinl 1 19991 ; iHerrnstein et al] 
l2005h . we can also establish a lower limit on s. Substituting 
the upper limit on Eo determined above into equation (9), 
we found that only for s > —2.3 the upper limit on the the 
disc mass is satisfied. As an example, we present in Figure 
1 the cumulative mass of the accretion disc as a function 
of radius, considering s — —2, —1 and 0. Note that this 
quantity resides below the observational upper limit for the 
disc mass of NGC 4258 in all three cases. 

3.2.3 Accretion rate, viscosity and bolometric luminosity 

The bolometric luminosity of a source is related to the mass 
accretion rate through: 



L bol = eMc 2 , (19) 
where 

e = 1 _ 1 ~ 2A ™ + a * A ™ 2 (20) 
\/l-3A m l + 2a«A^ 3 s /2 

is the accretion efficiency onto a Kerr black hole (e.g., 
IShapiro fc Teukolsky|[l983l) . 

The bolometric luminosity of NGC 4258 is roughly 

known, ranging approximately from 10 41 to 10 43 erg s _1 , 
which translates to 

1.765 x 1CT 6 . ■ ,„ _i, . !-765 x 1(T 4 

< M (M s yr l ) < . (21) 

However, there is a more restrictive upper limit for M 
since the condition £™ ln < E™ ax must be always satisfied. 
This relation leads to 
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M < 1.60 x l(T 4 a M Q yr" 1 , 



(22) 



wh ich is compatible with the r ecent estimate s obta ined 
by IHerrnstein et all i|2005h and iModiaz et all (120051 ). as 
well as with the a ccreti on rate inferred previously by 
iNeufeld fc Malonevl l|l995h . It is important to emphasize 
that this limit does not depend on the value of the power-law 
index s. 

On the other hand, the upper limit of M obtained 
from equation (22) must not be lower than its mini- 
mum value calculated from equation (21), implying that 
eat > 0.011, which translates to a lower limit for the 
a parameter (a > 0.026 for e < 0.42). Numerical 
magneto- hydrodynamic simulations, as well as observa- 
tions of many types of accretio n disc systems have in- 
dicated that 0.01 < a < 0. 2 JSiemiginowska fc Czernv 



1989; [Brandenburg et alj 119951; lHawlev. Gammie fc Balbus 



1996 



Jstl 



one et aU 119961; ^mak 1999; iMenou et all 



Hawlev fc Krolikl l200ll; IWinters. Balbus fc Hawlevl 1200 



and see King. Pringle fc Livid 2007 for a discussion) . 

Regarding the black hole spin, the conditions a < 0.2 
and ea > 0.011 are fulfilled only if a, > —0.07 in equation 
(20), which might be suggesting that the supermassive black 
hole in NGC4258 is rotating progradely in relation to its 
accretion disc unless it is spinning very slowly. 

In Figure 2 we show the parameters £0, tM, a and 
iboi/^Edd, where LEdd is the Eddington luminosity, as a 
function of the black hole spin for s —-2, -1 and 0. 

In Table 1 we summarize the physical parameters of 
NGC 4258 that will be used to analyze the feasibility of the 
Bardeen-Petterson effect in this particular source. 



4 BARDEEN-PETTERSON EFFECT IN 
NGC 4258 

4.1 Bardeen-Petterson radius 

Considering the extreme values of £0, ol and M for s = 
—2, —1 and (see Figure 2), we found in all cases that 
Rt 3> 10 7 -R ma , much greater than the outer radius of the 
maser disc. Therefore, if our disc parameterization is cor- 
rect, the accretion disc of NGC 4258 responds diffusively to 
the Bardeen-Petterson effect, which implies that we must 
use equation (3) to obtain _Rbp • As we mentioned in section 
2, we are assuming that U2 ~ f(a)ui; plugging it into equa- 
tion (3), we were therefore able to determine the Bardeen- 
Petterson radius for each set of disc parameters. 

In Figure 3 we show the Bardeen-Petterson radius as 
a function of the black hole spin for NGC 4258 for three 
different power-law surface mass density distribution (s = 
—2,-1 and 0). We have also considered four distinct cases 
related to different combinations among the parameters £0, 
a and M in order to produce the extreme values for Rbp- 

We can see that the values of the Bardeen-Petterson 
radius are comparable to or smalle r than 0.14 pc, the inne r 
radius of the maser disc inferred bv lHerrnstein et alj (|2005h . 
At first glance, this result may be seen as arguing against 
the Bardeen-Petterson effect since the warp has been in- 
ferred observationally in the outer disc. However, numerical 



simulations of Bardeen-Petterson discs, for several distinct 
physical conditions, usually indicate a smooth transition be- 
tween the inner (fiat) Bardeen-Petterson disc and the outer 
( misali gned) disc, exhibiting a g r adient in the inclination 
angle dNelson fc Papaloizoul 1 120001 : iFragile fc Anninosll2005l ; 
lLodato fc Pringlg|2006h . Consequently, warping in the outer 
disc could be associated with the gradient in the disc tilt an- 
gle, produced by the transition between a Bardeen-Petterson 
disc, and a misaligned outer disc that includes the maser re- 
gion. 

Indeed, the best-fitting model for the maser kinemat- 
ics is a warped disc with a radial inclination gradient 
l|Herrnstein et all l2005h . This model leads to a difference 
in inclination of about eight degrees between the inner and 
outer maser radius. Considering the full range of values for 
the Bardeen-Petterson radius displayed in Figure 3, and ex- 
trapolating the inclination gradient to the inner disc, we 
have found a variation in the tilt angle at the Bardeen- 
Petterson radius and innermost part of the maser disc that 
ranges from four to eight degrees approximately. 

Besides the warp in the tilt angle, there is also strong 
evidence that the maser disc of NGC 4258 is twisted az- 
imuthally, in the sense that the position angle in the 
disc varies with the radial distance from the black hole 
(e.g.. IHerrnstein et alj 20051 and references therein). Indeed, 
IHerrnstein et alj (|2005h found that a variation of about 10° 
in the warp position angle is sufficient to accommodate the 
observed maser kinematics. Assuming that the warping in 
the position angle is described by the quadratic function 
given by these authors, we find a variation of ~ 17° in this 
quantity from the inner maser radius to the extrapolated 
value at r ~ mas. Interestingly, a disc twisting of about 
10° is obtained if we consider the variation of the position 
angle from the core to the highest valu e of the Bardeen - 
Petterson radius calculated in this work. ICecil et all (|2000h 
and IWilson. Yang fc Cecil (|200ll ) proposed that jet preces- 
sion by an angle of about 10° can also reproduce the physi- 
cal characteristics of the anomalous arms. The similarity be- 
tween the required precession and warp position angle seems 
to suggest that they could be arising from the same physical 
mechanism. As we shall discuss in Section 4.3, the Bardeen- 
Petterson effect can provide the required jet precession in 
NGC 4258. We therefore propose the same mechanism to be 
responsible for the azimuthal twist of the disc. 

The allowed range for Rbp decreases for smaller values 
of the black hole spin, maintaining this trend even for retro- 
grade cases (as a consequence of the more restrictive limits 
on the disc parameters in those cases). 

Comparing the solutions represented by the full trian- 
gles and the full circles in Figure 3, we can infer the influence 
of a on Rbp ■ Fixing s and a, , the Bardeen-Petterson radius 
is larger when a is at maximum, increasing the coupling 
between consecutive disc annuli and consequently the effi- 
ciency of the Bardeen-Petterson effect. In fact, we can see 
from equation (3) that Rbp oc f(a) 1 ^ 3 ~ 1 \ for 7?bp 2S> Rg- 

The influence of £0 on _Rbp can be examined from a 
comparison between the open circles and full triangles: the 
Bardeen-Petterson radius increases with £0, which reflects 
the fact that _Rbp oc £^" 1,/ ' s-1 \ A comparison between the 
solutions represented by the full and open triangles shows 
that a decrease in the accretion rate results in an expansion 
of the Bardeen-Petterson radius, since 7?bp oc M 1 '' 8-1 ^. 
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Figure 3. The Bardcen-Petterson radius as a function of the 
black hole spin. From the upper to the lower panel, s varies from 
-2 to 0. Open circles correspond to the case of minimum value 
for So and M, and a = 0.2 (the maximum value considered in 
this work). Full circles represent the maximum value for So and 
the minimum for M and a. Full triangles refer to the case of 
minimum value for M and maximum for So and a, while open 
triangles represent the maximum for the three disc parameters. 
The dashed lines mark the position of the inner radius of the 
maser disc (r = 0.14 pc). 



Figure 4. The alignment time-scale as a function of the black 
hole spin for NGC 4258. From the upper to the lower panel, s 
varies from -2 to 0. As in Figure 3, open circles correspond to the 
case of minimum value for So and M, and a = 0.2, the maximum 
value considered in this work. Full circles represent the maximum 
value for So and the minimum for M and a. Full triangles refer 
to the case of minimum value for M and maximum for So and 
a, while open triangles represent the maximum for the three disc 
parameters. 



4.2 Alignment time-scale due to the 
Bardeen-Petterson effect 

In Figure 4 we present the alignment time-scale as a func- 
tion of the black hole spin, calculated from equation (5). 
For a fixed set of disc parameters, we can see a clear anti- 
correlation between T a i ign and s (in the sense that as s in- 
creases, the alignment time-scale becomes shorter), which is 
a consequence of the inertia of the disc (the disc mass inside 
the Bardeen-Petterson radius decreases as the value of s in- 
creases, as shown in Figure 1). The influence of a, So and 



M on the alignment timescale can be analysed (in a similar 
fashion to our discussion of the Bardeen-Petterson radius) 
by comparing the different symbols displayed in Figure 4 
(T a i ign oc So 1 ^ 1 / 2 -* , for Rbp > Jig). 

Our results show that the Bardeen-Petterson effect act- 
ing upon the accretion disc of NGC 4258 leads to alignment 
time-scales which are typically in the range from 0.4 to 10 
Gyr. 
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4.3 Anomalous arms and the radio jet of 
NGC 4258 

The anomalous arms, observed from radio to X-ra 

l ll 1 

wavelengths (e.g.. iBurbidge. B urbidgc & Prcndcrgastl 1963 

van der Kruit. Port. & Mathcwson 1972; ICecil et al. 200C 



Wilson. Yang fc Cecill 120011 ) , are stron gly believed to lie in 



the galactic disc of NGC 4258 (e .g. . Ivan der Kruitl fl974 



van Albada fc van der Hulsfl Il982l ; I Wilson. Yang fc Cecill 
20011 ). even though their nature has been rel ated to the 



jet p r opagation through the disc in some cases ([Ford et all 
I l98fj ; iMartin et all 1 19891 ; ICecil. Wilson fc Tullvlll992l ). The 
anomalous arms are straight at the inner parts of the disc 
and show, after about 2 kpc, a gradual curvature in oppo- 
site directions in relation to the nuclear region, forming a 
S-s hape structure. 

IWilson. Yang fc Cecill <|200ll ) performed a detailed 
study of NGC 4258, using high-resolution X-ray observa- 
tions, correlating these data with radio and optical images 
obtained in previous works. They proposed that the anoma- 
lous arms are shocked regions produced by the interaction 
between the disc gas and the mass motions induced by the 
propagation of the radio jet/counterjet. In addition, they 
suggested two different scenarios for the curva ture of the 
anoma lous arms: jet precession, also proposed bv lCecil et alj 
(2000), or buoyancy of the less dense and less tightly bound 
gas in the outer galactic disc. 

Here we shall assume that the curvature of the outer 
anomalous arms of NGC 4258 is driven by jet precession. 
From this assumption, we shall analyze the possibility that 
this precession is being produced by the Bardeen-Petterson 
effect. 

The Bardeen-Petterson effect induces not only a 
disc alignment but also a preces sion around the rota- 
tion axis of the black hole (e.g., Scheuer fc Feiler 1 19961 ; 



ICaproni, Abraham fc Mosquera Cuestal I2004T 120061 ). In- 
deed, Scheuer fc Feilerl 1 19961 ) showed that the alignment 
and precession time-scales are identical, with the precession 
period varying exponentially with time: 



P plec {At) = Tj 



align*- 



-At/T a i ign 



(23) 



where P prcc is the precession period after an elapsed time 
interval At. 

The precession angle Ai9 as a function of At and T a ii g n 
can be written as: 



A#(At) = 2tt 



At 



At/T align 



align 



(24) 



In order to be able to reproduce the physical charac- 
teristics of the anomalous arms, as well as the radio jet, 
it is necessary fo r the jet direction to chang e by ~ 10° 
l|Cecil et all l200d : IWilson. Yang fc Cecill l200ll ). which im- 
plies that A-& » 10° . Therefore, we can use equation (24) to 
estimate the time interval required to obtain a precession of 
the jet by ten degrees Atiodeg, for each solution derived in 
the previous sections. The results are shown in Figure 5. 

We note that a time interval between 5 and 300 Myr is 
necessary for the jet to change its direction by ten degrees. 
Such time-scales seem to be feasible since they are similar to 
or shorter than the viscous time-scale at 0.14 pc (between 
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Figure 5. Time interval for the jet/accretion disc of NGC 4258 
to precess by an angle of 10 degrees due to the Bardeen-Petterson 
effect as a function of the black hole spin. Symbols have the same 
meaning of those plotted in Figure 4. 



0.06 and 3 Gyr considering all combinations of the disc pa- 
rameters used here). Therefore, the Bardeen-Petterson ef- 
fect is an excellent candidate for driving the precession in 
NGC 4258. 



5 CONCLUSIONS 

In this work, we have analysed the viability of the Bardeen- 
Petterson effect as the mechanism responsible for the warp- 
ing and precession of the disc of NGC 4258. Evidence for a 
warped accretion disc in this source comes from the kine- 
matics of t he water maser infe r red from high-resolutio n ob- 
servations (jMivoshi et al.lll995l : lrTerrnstein et al.ll2005l ). Pre- 
cession of the radio jet/counterjet has been pointed as the 
potential reason for the curvature of the outer anomalous 
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arms in the galactic disc of NGC4258 (|Cecil et al.l l200d ; 
IWilson. Yang fc CecilllioOlh . 

Assuming a power-law surface density distribution for 
the accretion disc material, we used constraints inferred from 
observations and determined the acceptable ranges for the 
values of the accretion rate and the surface density at the 
inner radius of the disc. Our results also show a dependency 
of the black hole spin parameter on the power-law index 
of the surface density, such that a» cannot be smaller than 
about -0.07 if the observational constraints are to be obeyed. 

We have calculated the Bardeen-Petterson radius for 
a variety of disc parameters and, in all cases, it is similar 
to or located inside of the inner radius of the maser disc. 
This argues for the existence of a Bardeen-Petterson disc 
in NGC 4258. Consequently, the observed warping may be 
interpreted as resulting from a gradient in the disc inclina- 
tion angle due to the transition from an inner (flat) disc 
and a misaligned outer disc. All solutions lead to a time- 
scale for the alignment between the accretion disc within 
the Bardeen-Petterson radius and the black hole's equator 
shorter than about 10 Gyr. The alignment time-scales ob- 
tained were generally longer than about 0.4 Gyr. 

Assuming that the time evolu tion of the precession pe - 
riod follows the results obtained bv lScheuer fc Feilerl l|l996h . 
we calculated the time interval necessary for the accretion 
disc, and consequently for the jet/counterjet, to precess by 
an angle of ten degrees. This value is necessary to reproduce 
the curvature of the outer anomalous arms in the kiloparsec 
galactic disc of NGC 4258 if jet precession is responsible for 
it. We found that an interval between 5 and 300 Myr can 
make the jet to change its orientation by the required angle. 
This time-scale is shorter than the radial viscous time-scale 
calculated at 0.14 pc (between 0.06 and 3 Gyr, depending 
on the disc parameters), the inner radius of the maser disc. 
Given the similarity between the value of the jet precession 
angle and that of the twist of the azimuthal position angle 
inside the Bardeen-Petterson region of the disc, we suggest 
the possibility that both effects are associated with the same 
physical origin (frame-dragging by a Kerr black hole). 

We conclude that the Bardeen-Petterson effect may ex- 
plain the observations of NGC 4258 on both the sub-parsec 
and kiloparsec scales. 
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